Diffusion conductivity and weak localization in two-dimensional structures with 
electrostatically induced random antidots array 
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Results of experimental study of the weak localization phenomenon in 2D system with artificial 
inhomogeneity of potential relief are presented. It is shown that the shape of the magnetoconduc- 
tivity curve is determined by the statistics of closed paths. The area distribution function of closed 
paths has been obtained using the Fourier transformation of the magnetoconductivity curves taken 
at different temperatures. The experimental results are found in a qualitative agreement with the 
results of computer simulation. 



The transport properties of the semiconductor two- 
dimensional (2D) structures with antidots array were 
intensively studied both experimentally^'^'^'^'^'^'^'^ '^^i ^ ^ 
and theoreticallyi^ii^ii^i^ in last 10-15 years. The 
main concern was with the investigations of the baUistic 
systems, in which I > d, D, where / is the mean free path, 
d and D are the size of antidots and period of the antidots 
array, respectively. The rich diversity of the transport 
phenomena, such as commensurable oscillations, the pe- 
culiarities due to the trajectories roUing along the array 
of antidots was observed. The antidots array structures 
in the diffusion regime {I < d, D) were not studied essen- 
tially. Such structures are interesting in some aspects. 
First of all, the quantum corrections to the conductivity 
due to weak localization (WL) and electron-electron (e- 
e) interaction have to reveal the specific features when 
the phase breaking length, — ^JVt^, where V is the 
diffusion coefhcient and is phase breaking time, or the 
temperature length, Lt = yP/T, become larger than 
d, D at decreasing temperature (hereafter we set fc^ = 1, 
h — 1). Secondly, the large enough negative gate volt- 
age has to deplete the channels between the antidots (the 
channel width wq is about D—d) and as a result to lead to 
crossover to the hopping conductivity. Besides, from the 
transport properties standpoint the antidots arrays are 
the fine model of the granular media. In contrast to the 
granular metallic film, the parameters of the "granules" 
and "barriers" are reliably known and can be changed 
continuously within wide range. 

In this paper we report the results of the experimen- 
tal study of the weak localization correction to the con- 
ductivity in the structure with the random array of an- 
tidots. We show that the change of the magnetoresis- 
tance at arising of the antidots and increase of the an- 
tidots size results from the change of statistics of closed 
paths. Namely, the contribution of the trajectories with 
the large enclosed area is strongly suppressed. The ex- 
perimental area distribution function is in a reasonable 
agreement with the function obtained from the computer 
simulation. 

The random antidots array was made on the single 
quantum well heterostructure with the electron density 
n = 1.5 X 10^2 cnr-2 and mobility = 19000 cmV(Vs) 
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FIG. 1: The sketch (a) and electron microscope image (b) of 
the structure investigated. 



grown by the molecular beam epitaxy. It consists of a 
250 nm-thick undoped GaAs buffer layer grown on semi- 
insulator GaAs, a 50 nm Alo.3Gao.7As barrier. Si 5 layer, 
a 6 nm spacer of undoped Alo.3Gao.7As, a 8 nm GaAs 
well, a 6 nm spacer of undoped Alo.3Gao.7As, a Si J 
layer, a 50 nm Alo.3Gao.7As barrier, and 150 nm cap layer 
of undoped GaAs [Fig. [IJa)]. The samples were etched 
into standard Hall bars. The holes in the cap layer [see 
Fig. Hl^b)] were fabricated with the use of electron beam 
lithography and wet etching. Their depth measured by 
atomic force microscope (AFM) consists of about 85 nm. 
Note that this depth is less than the cap layer thick- 
ness before etching. The holes of 0.7-/xm-diameter were 
shifted randomly on the value ~ 0.1 jum from the sites of 
the square lattice with the period of about 1 fim. This 
shift destroys all oscillations in the galvanomagnetic ef- 
fects resulted from the commensurability between the cy- 
clotron orbits and the lattice period and from the trajec- 
tories rolling along the array of antidots. After etching an 
Al gate electrode was deposited by thermal evaporation 
onto the cap layer. 

The magnetic field dependences of the longitudinal re- 
sistance {pxx) for the pattern structure together with that 
for the unpattern one for some gate voltages are shown 
in Fig. [21 It is seen that they are rather complicated. 
Following the sharp magnetoresistivity in low magnetic 
field, which results from the suppression of the interfer- 
ence quantum correction, the relatively smooth negative 
or positive magnetoresistivity against the background of 
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FIG. 2: The magnetic field dependences of pxx for pattern 
(left column) and unpattern (right column) structures for T = 
4.2 K. The dotted lines in (b) and (c) are calculated from 
the standard classical formula for two types of carriers with 
the electron densities found from Eqs. ^ and ([2}, and /ii = 
12600 cmV(Vs), /ia = 3400 cmV(Vs) for Vg = -1 V, and 
/ii = 8600 cmV(Vs), = 1300 cm7(Vs) for Vg = -2.4 V 
that corresponds to the model described in the text. 



the Shubnikov-de Haas (SdH) oscillations is observed. 

Let us first consider the range of high magnetic field. 
At Vg = the magnetoresistance in the pattern and 
unpattern structures are very close to each other [see 
Figs.mja) and[2{f)]. In both cases the paraboliclike nega- 
tive magnetoresistance resulting from the e-e interaction 
correction (for more detail see Ref. [13) and the SdH oscil- 
lations of close frequency are observed. Such the behavior 
is not surprising since the difference in local electron den- 
sities in quantum well under the holes and out of them, 
6n, is relatively small at Vg = 0. It can be estimated as 
5n ~ (C2 — Ci)Vs/\e\, where Vs is the surface potential, 
Ci,2 = ££0/^1,2 is the local capacity, ti,2 is the distance 
between 2D gas and gate electrode in different locations 
of the structure [see Fig. [Ha)], eo is dielectric constant 
of free space. With the use of Vs — 0.7 V, ti = 85 nm, 
^2 = 125 nm, and e = 12.5 we obtain 5n/n w 0.1. How- 



ever, 5n strongly increases with the lowering gate voltage 
that leads to the positive magnetoresistance evident in 
the pattern sample within the gate voltage range from 
— 1 V to —1.8 V. The rise of the positive magnetoresis- 
tance is a sequence of the comparable contributions to 
the conductivity from the regions under the holes and 
out of them. For the first approximation the transport 
in such inhomogeneous media can be considered as deter- 
mined by two types of carriers with the different mobility 
and densityii^ The dotted curves in Figs. [2{b) and O^c) 
are the magnetoresistivity calculated in the framework 
of this simple model with parameters determined below. 
At Vg < —1.8 V, the negative magnetoresistance is re- 
stored because the dielectric holes in 2D gas are formed 
and the conductivity of the structure is determined by 
the channels between the antidots. 

As figure [2] shows the SdH oscillations in the pattern 
structure are observed down to Vg = — 2.8 V. The peri- 
ods of the oscillations in the pattern and unpattern struc- 
tures at given Vg are close to each other [see Fig. (Ha)] 
and the gate voltage dependence of the electron density 
calculated from the oscillations is well described by the 
expression 

ni{Vg) = n{Vg) = (1.52-1-0.33^9) X 10l^ cm^^ (1) 

The facts that the oscillations of only one period are ob- 
served in the pattern sample and the dependences ni{Vg) 
and n{Vg) are practically the same mean that only the 
areas of 2D gas located out of the holes contribute to 
the SdH oscillations. The areas under various holes have 
probably different electron density due to the different 
depths of the holes and, therefore, the corresponding os- 
cillations are very broadened. The Vg dependence of the 
electron density under the holes n2{Vg) can be obtained 
from the geometric consideration. With the use of the 
local cap layer thickness t2 — 125 nm we have 

n2{Vg) = (1.35 + 0.55 Vg) x 10^^ cm'^. (2) 

In Fig. IHa) this dependence is shown by dashed line. 

Thus analyzing the high-field magnetoresistance we 
reason that: (i) at < —1.8 V the conductivity is 
mainly determined by the channels; (ii) the electron den- 
sity out of the dielectric holes remains more or less ho- 
mogeneous in this Vg range. 

Let us compare the behavior of the conductivity for 
pattern {opatt) and unpattern {aunpatt) samples at i? = 
0. The values of apatt and (Jimpatt measured at T = 4.2 K 
are plotted against the electron density found from the 
SdH oscillations in FiglJfb). It is seen that the conduc- 
tivity of the pattern structure significantly steeper falls 
down with decreasing n than that of unpattern sample. 
Qualitatively such a behavior is transparent. This is be- 
cause that the 2D gas under the holes in the pattern sam- 
ple is depleted faster with lowering Vg than that out of 
them due to thinner cap layer in these locations. If this is 
the case we can obtain the geometrical parameters of the 
conducting areas knowing the experimental value of the 



3 



V(V) 

-3 -2 -1 






1 ' 

LlrL 


1 ' 1 ' 1 ' 1 ' 
























y * 






/ 






p ' 






/ ' 






r * 






/ < 






/ p 






/ f 






5 r 






1 






n 






( 

J.I.I 



0.5 1.0 1.5 
n (10''cm ') 



FIG. 3; (a) The gate voltage dependence of the electron den- 
sity found from the SdH oscillations for the pattern (solid 
symbols) and unpattern (open symbols) structure. Solid and 
dashed lines are drawn according to Eqs. ([1} and ((2)1, respec- 
tively, (b) The electron density dependences of the conduc- 
tivity measured at T = 4.2 K for the pattern (solid symbols) 
and unpattern (open symbols). 



ratio a unpatt I <y patt referred below as K. Using the data 
obtained for the unpattern structure we obtain the mean 
free path, I ~ 0.03 . . . 0.3 /zm depending on the gate volt- 
age, being less than the characteristic scales of the holes 
and channels, ^ 0.5 /xm [see Fig. [Ijb)]. Therefore, one 
can deal with the local conductivity. If one additionally 
neglects the randomness in the antidots position, we can 
write out the following approximate expression for the 
conductivity of the pattern structure: 



Cpatt 




dy 



criw{y) + cr2[D - w{y)] 



(3) 



where ai and (T2 stand for the local conductivity of the 2D 
gas out of and under the holes (in the insert of Fig. [D^a) 
these areas labeled as 1 and 2, respectively), w{y) is the 
y-dependent width of the area 1. In what follows we 
suppose the area 2 being round in the shape. The above 
equation gives the result, which coincides with the exact 
solution with the accuracy better than 20 percent when 
wq/D > 0.1, where wq = w{D/2). In order to calculate 
the dependence KiVg) it is natural to suppose that the 
local conductivity CTi and a2 is fully determined by the 
local electron density ni and n2, respectively, and the 
ai-vs-Ui behavior is just the same as that for the unpat- 
tern sample cri(ni), CT2(n-2) = (^unpatt{n) [shown by open 
symbols in Fig. [3l[b)]. 

In Fig. mja) we present the K-ys-ni dependences as 
they have been obtained experimentally and calculated 
from Eq. ^ with d = 0.7 /im obtained from AFM. It 
is seen that the above simple model well describes the 
experimental results down to rii ~ 8.5 x 10^^ cm~^ that 
corresponds to ~ —2 V. The reason for the discrep- 
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FIG. 4: (a) The K-YS-ni dependence obtained experimentally 
(symbols) and calculated from Eq. ([3]) (lines) for d = 0.7 /xm 
and D — 1 /im. Insert shows one period of the antidots ar- 
ray, (b) The ni dependence of the channel width wq obtained 
from the experimental values of K at low Vg with the help of 
Eq. (131). The line is provided as a guide for the eye. 



ancy which is clearly evident at lower electron density is 
transparent. For the fixed d value, the saturation of the 
calculated -fiT-vs-ni dependence at ni < 8.5 x 10^^ cm~^ 
{Vg < —2 V) results from the fact that a2 becomes much 
less than ai. The enhance of K obtained experimen- 
tally for these gate voltages is sequence of the depletion 
of the area outside the antidots, i.e., of an increase of 
the antidots size d and decrease of the channel width 
wq = D — d. Thus knowing the experimental value of K 
and using Eq. Q we are able to find the channel width 
Wq when ni < 8.5 x 10^^ cm~^ . The results are depicted 
in Fig. lUb). It is seen that the separation between an- 
tidots Wo decreases with ni decrease. Extrapolating the 
WQ-vs-ni plot to Wo = one obtains that the antidots 
close when ni ~ 4 x 10" cni^^ [Vg = -(3.4 . . . 3.2) V]. 

Thus, the gate voltage dependence of the conductivity 
and the high magnetic field magnetoresistance are rea- 
sonably described within the following simple model. At 



~ (0 ... — 2) V the conductivity is determined both by 
the areas under holes and out of the them. These areas 
are characterized by the different electron density, which 



determines the local conductivity. At Vg 



.3) V 



the antidots are formed and the conductivity of struc- 
ture is determined by the channels between the antidots 
with local conductivity cti equal to the conductivity a of 
the unpattern sample at the same electron density. Fi- 
nally, at Vg = —(3.4 . . . 3.2) V the channels are collapsed 
and most likely the crossover to the hopping conductivity 
should occur. 

Let us now inspect the low magnetic field negative 
magnetoresistivity which results from suppression of the 
WL correction. We focus our consideration on the re- 
sults obtained within the second range of the gate volt- 



ages: Vg 



-(2... 3) V. The magnetic field depen- 
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FIG. 5: (Color online) The low-field magnetoconductivity for 
different temperatures for the pattern (a) and unpattern (b) 
structure. The local conductivity is shown in panel (a). 



dences of the local conductivity measured in units of 
Go = (?I{2'k'^K) ~ 1.23 X 10^3 Q-i are presented for 
different temperatures in Fig. [51[a). For comparison, the 
analogous dependences for the unpattern structure mea- 
sured at close conductivity value are shown in Fig. [SJb). 
For the first sight the magnetoresistance curves in the 
panels are very similar. However, this impression is 
wrong. The difference in magnetoresistivity shape for 
these structures is more pronounced when comparing the 
results of data treatment performed in a standard man- 
ner. 

The shape of low-field positive magnetoconductivity 
A(t{B) = p~^{B) — ^^2^(0) caused by suppression of the 
weak localization in homogeneous 2D gas is described by 
the Hikami-Larkin-Nagaoka (HLN) expressioni2i2£ 



A(t(B) 



B 



\Ttj, Btr , 



lux, (4) 



where Btr = ft/ {2eP) is the transport magnetic field, r 
is the momentum relaxation time, tpi^) is a digamma 
function, and a is the prefactor, which is equal to unity 
in the diffusion regime {B < Btr, t t^) and at high 
conductivity (ct ^ Go). 

We have used this expression to fit each experimental 
curve within different range of magnetic field with a and 
T0 as the fitting parameters. The results are shown in 
Fig. [SI How the fitting parameters and a depend on 
the range of magnetic field is shown in Figs. [6ja) and 
ini^b), whereas their temperature dependence is shown in 
Figs, ^c) SLiid^d). One can see that the parameters 
found for the unpattern structure behave themselves rea- 
sonable. Their values only slightly depend on the fitting 
interval. The prefactor is close to unity and practically 
independent of temperature. The temperature depen- 
dence of is close to the theoretical one oc 1/T. 
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FIG. 6; The magnetic field and temperature dependences of 
the fitting parameters and a for pattern (solid symbols) 
and unpattern (open symbols) structures. The data in (a) 
and (b) are shown for the lowest temperatures, the data in 
(c) and (d) are obtained for the fitting interval from up to 
0.25Btr. 



Thus, the HLN expression well describes A(j{B) for the 
unpattern sample. 

In contrast to this, the strong sensitivity of the fitting 
parameters to the fitting interval takes place for the pat- 
tern sample. Therewith the value of a is significantly 
larger than unity and strongly dependent on the tem- 
perature. The value of the fitting parameter is much 
less than that for the unpattern structure and it satu- 
rates with the decreasing temperature. All of this means 
that the WL correction for the pattern structure is not 
described by Eq. (j4]) and the determination of the phase 
breaking time by the standard way is impossible.— 

The strong dependence of the fitting parameters on 
the fitting interval of the magnetic field indicates that 
the role of dielectric antidots in 2D gas is not reduced 
to arising of the prefactor exp(— t^/t^) with te as the 
Ehrenfest time, which suppresses the weak localization 
in the 2D gas with hard discs as scatterer o^'^^i^^ . 

Specific features of the weak localization in the pattern 
sample can be understood by considering the quasiclassic 
interpretation of this phenomenon. Within quasiclassic 
approximation the conductivity correction is expressed 
through the classical quasiprobability for an electron to 
return to the area of the order Xpl (Xp = 2Tr/kF, kp is 
the Fermi wave vector) around the start poinij2^i^i21i2S 



5(7 



(5) 



where cto = T^kplGo is the Drude conductivity, and W 
stands for the quasiprobability density of return {quasi- 
means that W includes not only the classical probability 
density, but the interference destruction due to an ex- 
ternal magnetic field and inelastic scattering processes). 
With taking these effects into account Eq. ^ can be 
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FIG. 7: (Color online) The area distribution function W{S) 
obtained experimentally for the unpattern samples for two 
gate voltages (symbols). Line is the diffusion limit: Anl^W = 
1/5. 



rewritten as follows2S 

Sa{b) = -27rZ2Go / dS \ w{S) 

J — OO L 

m 



exp 



(6) 



where is the phase breaking length connected with 
through the Fermi velocity, — vpr^, W{S) and 
L{S) are the algebraic area distribution function of closed 
paths and the area dependence of the average length of 
closed paths respectively (for more detail see Sec. II of 
Ref. i2&) . This equation shows that the shape of the mag- 
netoconductivity curve is determined by the statistics of 
the closed paths, namely by the area distribution func- 
tion W{S) and by function L{S). It is clear that the exis- 
tence of dielectric dots in 2D gas should change the statis- 
tics of closed paths resulting in the change of the shape 
of the magnetoconductivity curve. In Ref. [sol . there is 
shown how the analysis of the Fourier transform of the 
negative magnetoresistance provides the information on 
the function W{S). The short of the matter is clear from 
Eq. ([6]). It is seen that the Fourier transform of 5a{B) 



1 

$(5*) = / dB 6a{B) cos 



is equal to 



$(5) = -2TTl^GoW{S)e^p 



f2TTBS\ 



L{S) \ 
h ) 



(7) 



where $0 = 7rS/e is the elementary flux quantum. Since 
tends to infinity when T — > 0, the extrapolation of 
$(S', T) to T = should give the value of 27r;2Gol^(S'). It 
is ideal situation. In the reality, such the approach allows 
us to obtain experimentally the area distribution function 



FIG. 8: (Color online) (a) The function ATrt^WiS) obtained 
experimentally for the pattern sample for different gate volt- 
ages. Arrows indicates where the plots deviate from 1/S- 
dependence. (b) The function A-kI'^W{S) obtained from the 
simulation procedure carried out with parameters correspond- 
ing to the gate voltages in panel (a). 



of the trajectories which length is less than approximately 
21'^, where l'^ is the phase relaxation length at lowest 
temperature of the experiment, because the contribution 
of longer trajectories to the magnetoconductivity is very 
small. 

The results of the data processing for the unpattern 
and pattern samples are shown in Figs. [7]and[Hfa), re- 
spectively (note the scales in these figures are identi- 
cal). It is seen that not only the W^-vs-S" dependences 
are drastically different for the pattern and unpattern 
samples but their responses to the change of the gate 
voltage as well. The area distribution function for the 
unpattern sample is very close to that predicted theoreti- 
cally for 2D homogeneous systems in the diffusion regime, 
At:PW = S-^, S>P (see Refs. [li andllO), it is practi- 
cally insensitive to the gate voltage (Fig. [7]). As for the 
pattern sample, it is seen that the curves follow near the 
\/ S dependence only at low S < Sc values. At higher ar- 
eas, S > Sc, W decreases with S increase much steeper 
than S~^; the lower gate voltage is the lower areas at 
which this deviation occurs [Fig. [8l^a)]. 

Before to interpret such the behavior of W{S) for the 
pattern sample let us clarify the meaning of this quantity 
for this case. Because our system is inhomogeneous due 
to dielectric inclusions under the holes in the cap layer, 
Eq. ^ should be rederived. This is because the function 
W{S) is not universal now, it should depend on the po- 
sition of starting point. If one neglects as above the ape- 
riodicity of the antidots array, we obtains the following 
expression connecting the correction to the conductivity 
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of the antidots array and local correction Saf. 
f _dy 

Sa{b) = 



dx5ai{x, y, b) 



dy 



1 2 



wiy) 



(8) 



Here, the integration runs over the intervals given by the 
border of conducting area, wly) is the y dependence of 
width of conducting area, and the quantity ai{x,y,b) is 
given by the expression analogous to Eq. 



6(Ji{x,y,b) = -27r/2Go I dS { 

J —OO L 



mx,y,S) 



exp 



L{x,y,S) 

Id, 



cos 



(9) 



in which W{x, y, S) is determined in such a way that 
W{x, y, S)dS gives the density probability of return to 
the starting point with the coordinates {x, y) with the 
enclosed algebraic area in the interval {S, S + dS). 

Thus, the above described procedure, which for the 
homogeneous 2D gas allows us to obtain the area distri- 
bution function of closed paths, being applied to the data 
obtained for the pattern sample gives the effective area 
distribution function: 



W{S) 



dy 



[w{yW 



W(x, y, S) dx 



dy 
w{y) 



(10) 



Now we are in position to discuss the behavior of W{S) 
for the pattern sample [Fig. [51Ja)] . In order to understand 
main features we have performed the computer simula- 
tion of a particle motion over 2D plane with scatterers. 
The details can be found in Refs. [29ll3lll3^ . below is out- 
line only and important features. The 2D plane is repre- 
sented as a lattice with scatterers of two types placed in 
a part of lattice site. The scatterers of the first type with 
isotropic differential cross-section correspond to ionized 
impurity. The scatterers of the second type are hard discs 
with specular reflection from the boundaries. Particle 
motion is forbidden within the disks. They correspond 
to the areas of 2D gas under the holes. A particle is 
launched from some point with x, y as coordinates, then 
it moves with a constant velocity along straight lines, 
which happen to be terminated by collisions with the 
scatterers. After collision it changes the motion direc- 
tion. If the particle passes near the starting point at 
the distance less than some prescribed value a/2 ^ I, 
the path is perceived as being closed. Its length and en- 
closed algebraic area are calculated and kept in memory. 
The particle walks over the plane until the path traversed 
is longer than 21'^, where l'^ is the phase relaxation time 
obtained at lowest temperature on the unpattern sam- 
ple for the same Vg value. As mentioned above namely 
the statistics of such closed paths can be reliably ob- 
tained from the weak localization experiments. When 



the path becomes longer than this value another particle 
is launched and all is repeated. After large number of 
launches from the starting point with given x, y coordi- 
nates 2U(x,?/, S) is calculated as 



2U(x,y,5) = 



ns 



NlaAS' 



(11) 



where N is the number of starts from this point, ns is 
the number of returns along the trajectory with enclosed 
area in the interval {S, S + AS). Launching the particle 
from different starting points, we are able to calculate 
numerically the function W{S) using Eq. (fTO| in the dis- 
crete form. 

Shown in Fig. [Hljb) are the results of simulation pro- 
cedure carried out with the parameters corresponding to 
the gate voltages from Fig.[HIa). It is clearly seen that the 
results of computer simulation are in qualitative agree- 
ment with those of real experiment. 

Analysis of the simulation results gives an insight into 
why the experimental W-vs-S dependence is steeper for 
the pattern sample. This can be done if one consid- 
ers how the paths with the given area enclosed are dis- 
tributed over the length. This distribution is character- 
ized by the function ws{L) determined in such a way that 
ws{L)dS gives the density probability of return along 
a trajectory with the length L and area in the interval 
{S, S + dS). In Fig. [H] we show the functions ws{L) ob- 
tained from the simulation procedure for the pattern^ 
and unpattern samples (all the paths including those for 
which L > 21'^ are taken into account here) . Qualita- 
tively, the functions are similar. They have a peak which 
position characterizes the typical length of closed paths. 
In both cases, the larger area enclosed the longer trajec- 
tories. The ws-vs-L plot for the unpattern sample is well 
described by the expression 



wsiL) 



1 



[21 Lf 



■ cosh 



■kS 

Tl 



(12) 



obtained analytically within the diffusion approximation 
for the ordinary 2D system^^ that justifies the validity 
of simulation procedure. An important point is that the 
closed paths in the pattern sample are significantly longer 
than the paths in the unpattern sample with the same 
area enclosed. If, in correspondence with the experimen- 
tal situation, one restricts the consideration by short tra- 
jectories, L <2l',, we obtain 



WiS) = 1^'' ws{L)^- 



(13) 



approximately the same for the both samples when S is 
relatively small [Figs.[9l[a)], and much less in the pattern 
sample when S is sufficiently large [Fig. OJc)]. Thus, the 
different length distribution of closed paths is the reason 
of the different behavior of the area distribution function 
obtained experimentally for the pattern and unpattern 
samples. 
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L(^m) 

FIG. 9: (Color online) The length distribution functions 
of closed paths with different area enclosed obtained from 
the simulation procedure with parameters corresponding to 
the unpattern (circles) and pattern (squares) samples for 
Vg = —2.84 V. Lines are calculated from Eq. (|12|l with cor- 
responding S-values and I = 0.06 /im. Hatched areas are the 
W value for paths with L < 2lL. 



Let us return to the magnetoconductivity caused by 
suppression of the interference quantum correction. The 
simulation procedure allows us to calculate 5a{h) for the 
model system. This can be done with the use of Eq. ([8]) 
in the discrete form and of 5ai{x^ y, b) calculated as 

where summation runs over all closed trajectories among 
a total number of trajectories Ni starting from the point 
with Xi and yi as coordinates, 5^ and stand for alge- 
braic area and length respectively of the fc-th trajectory. 

In Fig. [To] we compare the simulation and experi- 
mental results. Symbols are the experimental plot ob- 
tained for the pattern sample at = —2.25 V. Solid 
lines are the simulated magnetic field dependences of 
Act = 5(j{B) - 5(j{Q). The values of = 1.9 x 10"" s 
for T = 1.35 K and = 6.8 x IQ-^^ s for T = 4.2 K 
correspond to the best accordance with the experimental 
data. It is significant that they are close to that obtained 
for the unpattern sample [see Fig. EUc)]. If one calcu- 
lates Acr(i?) from Eq. (HJ with the same parameters and 




0.00 0.02 0.04 0.00 0.02 0.04 



B(T) B(T) 

FIG. 10: (Color online) The low-field magnetoconductance 
for the pattern sample for T = 1.35 K (a) and 4.2 K(b), Vg = 
—2.25 V. Symbols are the measured dependences. Solid and 
dashed lines are the simulation results and Eq. (j4]| obtained 
with = 1.9 X 10"" s (a) and = 6.8 x 10"^^ s (b). 

a = 1 we obtain drastic disagreement with the experi- 
mental magnetoconductance (see dashed lines Fig. [T0|). 

Thus, the simulation approach allows us to understand 
qualitatively the main features of the low-field magne- 
toconductivity in the pattern samples, which are deter- 
mined by peculiarities of statistics of closed paths. 

It would serve no purpose to make more detailed com- 
parison between the experimental and simulation results 
because our model is rather crude. So, we supposed that 
the areas forbidden for classical motion had the form of 
hard disks. As seen from Fig. [IJb) it is not exactly. More- 
over, it was suggested that the mean free path and Fermi 
velocity are constant over all the conducting area that 
should not be fulfilled in the real structures. 

In summary, we have studied the weak localization 
in the 2D electron gas with potential electrostatic relief 
forming the insulating antidots array. It has been shown 
that the use of the standard procedure for the obtaining 
of the phase relaxation time from the shape of the mag- 
netoconductivity curve is inadequate when it is applied 
to the antidots array structure in the diffusion regime. 
To understand the main features of the weak localization 
in these systems the alternative approach based on the 
analysis of the statistics of closed paths has been used. 
We have shown that the main peculiarities of the WL 
phenomenon in the pattern structure are due to the spe- 
cific statistics of closed paths. The paths of actual areas 
in system with antidots are characterized by significantly 
larger lengths as compared with the usual 2D systems. 
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